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. ' The probability of simultaneous occurence of at least k span- 

j^ I ning clusters has been studied by Monte Carlo simulations on 

C^ ' the 2D square lattice at the bond percolation threshold Pc = 1/2. 

I I It is found that the probability of k and more Incipient Span- 

' O ' ning Clusters (ISC) has the values P{k > 1) « 0.00658(3) and 

"' P{k > 2) « 0.00000148(21) provided that the limit of these prob- 

abilities for infinite lattices exists. The probability P(k > 3) of 
more than three ISC could be estimated to be of the order of 10~^^ 
and is beyond the possibility to compute a such value by nowdays 
computers. So, it is impossible to check in simulations the Aizen- 
QQ ' man law for the probabilities when k » 1. We have detected a 

^^ I single sample with 4 ISC in a total number of about 10^*^ sani- 

es I pies investigated. The probability of single event is 1/10 for that 

number of samples. 
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It was a common belief until a very recent time that on 2D lattices at perco- 
lation threshold pc there exists exactly one percolation cluster |j|, |2| . Indeed, 
'^ I it was rigorously proven by Newman and Shulman ^ that the number of 

Infinite Clusters is either 0, 1 or oo and later is was proven ||^, ^ that the 
O ' infinite cluster is unique. 

New insight developed recently by Aizenman, who proved Q (see also, 
his talk at StatPhys 19 P]) that the number of Incipient Spanning Clus- 
ters (ISC) in 2D critical percolation can be larger than one, and that the 
probability of at least k separate clusters is bounded 



PLik){ - ''',. (1) 



< 



where a and a' are positive and L is a linear lattice size. Moreover, he 
conjectured the existence of the limit 



hm hm 7^ log Piik) = aasymp ■ 
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Indications of the existence of simultaneous clusters in two dimensional 
critical percolation in the limit of infinite lattices were found in computer 
simulations of Parongama Sen g, ^ for site percolation on square lattices [| 
with helical boundary conditions and in a strip geometry by Hu and Lin 



[10|. The detailed history of the recent development can be found in recent 
Stauffer mini-review []TT| . 

Here we report on our attempt at a direct verification of Aizenman rig- 
orous result by computer simulations. 

We investigate by Monte-Carlo the number of spanning clusters in the 
critical bond percolation model on two dimensional square lattices. We have 
determined the first limit in (2), i.e. the numerical values of the probabilities 

P{k) = lim Piik) (3) 

L— >oo 

for A: = 1, 2 and 3, as shown in Fig. 3. 

Using calculated probabilities we can conclude that the probability P{k = 
4) of four ISC is of the order of 10~^^ and is beyond the power of today's 
computers. So, the numerical check of Aizenman result is a good task for 
the computers of XXI century. 

In the Section 2 we define precisely the model and the algorithm used 
in simulations. The influence of subtle details of alghorithms on the values 
of spanning probabilities was emphasized recently by Aharony and Stauffer 
[12|. In Section 3 we present the details of simulations. The discussion of 



the results is in the last section. 



2 Model and Algorithm 

We use in simulations a rectangular square lattice with linear size L, and 
exactly L sites and L bonds both horizontally and vertically. We use free 

^Actually, p. Sen found more than one spanning cluster not only in dimension 2 but 
also in dimensions 3,4 and 5. 



boundary conditions. For clarity the example of lattice is shown in Fig. 1 
together with the dual lattice. The dual lattice has the same number of 
sites and bonds as the original lattice (compare with [O] and 0]). This 
gives the possibility to keep in the finite lattices some properties of infinite 
lattice. First of all, the number of bonds is exactly twice the number of 
sites. Second, the self-duality is valid for any finite lattice size. Third, the 
horizontal and vertical directions are equivalent. 

Each bond could be occupied with probability q = 1—p and closed with 
the probability p. Given the realization to each of 2-L^ bonds to be occupied 
formed the sample. Each sample could be decomposed in clusters of con- 
nected occupied bonds. For that we use Hoshen-Kopelman |1J] algorithm, 
which is exact. We are interested in the spanning properties of such a clus- 
ters. Namely, what is the probability that a cluster connects the opposite 
borders of square and what is the number of disjoint spanning clusters? 

An event h that the cluster spans the lattice horizontally is an event that 
at least one of the left sites and at least one of the right bonds are in the 
same cluster. The probability of such event is just the Langlands et all [^] 
horizontal crossing probability vr/j. 

For our purposes we need more detailed events. Namely, by hi we will 
denote the event that there are exactly one cluster connecting left sites with 
the right occupied bonds. In the same manner by /i^ we will denote the event 
that there are exactly k disjoint clusters connecting left and right borders of 
our lattice. Then, the horizontal crossing probility vr/j is given by X^i^i '^hi i 
where we denote by nh^ the probability of event h^ ■ 

In a full analogy, by vr^j, we will denote the crossing probabilities from top 
to bottom, i.e. vertically. Obviously, for our choice of the lattice, iTh^ = tTv,, 
for any lattice size. This could be used as a check of the calculated statistical 
errors. 

Knowing the origin of the bad properties of both main fast methods for 
random number generation [|l5|, ll^ we check the results using the same linear 
congruential method as Langlands, et all use in [^] (xj+i = (a Xj+c) mod m, 
with a = 142412240584757, c = 11, m = 2^^) and the shift register (x„ = 
a;n-i57 ffi Xn-314 ® Xn-471 © 2;„_9689) (H^l (the One uscd by Ziff in [|l8|) and 
found that the results coincide, which suggests the absence of systematic 
errors. 

It should be noted that the variance of probabilities is independent from 
tha lattice size because probabilities are calculated as the expectation val- 
ues of corresponding indicator functions. So, we should keep the number of 
samples independent of the lattice size and, therefore, the only parameter 



which controls the statistical errors is the number of samples M. Probabil- 
ities were calculated averaging M = 10^ samples and the error bars were 
defined from 100 bins, each bin being the average over 10^ samples. 

Throughout the next section we deal with the critical bond percolation 
on square lattice {p = Pc = 1/2) with free boundaries. 



The simulations was done under the Topos environment [21| working on 
a number of Digital Alpha workstations and servers in Landau Institute and 
Chernogolovka Science Park. 

3 Numerical results 

We calculate the probabilities -nh^ and 7r^,^. for square lattices with linear 
sizes L = 8, 12, 16, 20, 30, 32, 64 looking for ah events up to 8.0 Within the 
statistical errors the vertical and horizontal crossing probabilities coincinde 
giving additional confidence in the quality of data (Table I). 

On the Fig.l we plot the probability for simultaneous occurence of more 
than one spanning cluster Pl(/c > 1) = X]fc>i '^h^ versus the inverse system 
volume 1/L^. A linear fit gives us the limiting value oiP{k > 1) = 6.58- 10~^ 
with error ~ 3 • 10~^. 

Fig. 2 shows the dependence of the probability of more than two clusters 
versus 1/L^. The best linear fit gives the limiting value of P{k > 2) = 
1.48 • lO"*^ with uncertainty 2.1 • lO"'^. 

Actually, we observed in computations mostly the events of up to three 
simultaneous spanning clusters. We simulated at total about 10^'^ samples of 
different sizes 8 < L < 64 and only one sample with 4 spanning clusters was 
detected. This event clearly not contradicts with our estimate for -P(4) ~ 
10~^^, given below. Single event is probable as one part in ten. 

We could fit exponent a in (2) taking q^ = — 1/A;^ log Pin > k) and 
obtaining ai = 0.693, 02 = 1.256 and 03 = 1.498. 

The logarithm of P{k) is plotted on Fig. 3 versus k"^ together with the 
best linear fit, which gives a oc 1.61(7) and ^ oc 3 it 1.5. For comparison, 
a linear fit for the logarithm of P{k) versus k gives the value of exponent 
2.8(5) with the too large uncertainty. 

This is as full analysis as we could do because even with the using the 
of a = 1.5 the probability of four ISC would be 3.8 • 10"^^ and of five ISC of 
5.5-10^^^. Even the first value is impossible to check with the computers one 



^This does not imply that we seriously expect any events with k > A. The upper limit 
of eight is simply related to the numerical algorithm. 



have today. To estimate accurately the probabihty of an event expected to 
be of the order 10~^^ one needs to generate at least 10^^ samples. Suppose, 
we have a computer with the CPU cycle of A ns and an algorithm which 
needs m CPU cycles to process one lattice site. Thus, it is m • L^ cycles 
per sample and 10^^ ■ m ■ L^ cycles altogether. The time needed will be 
10^^ ■ m ■ L"^ ■ A. If we are very optimistic, we could think that m = 10 and 
A = 1 ns, and taking a moderate lattice size L = 32 one gets about 1157 ... 
days! 

In addition, we compute the distribution of cluster sizes. In Fig. 4 we 
plot the mean value < s > per lattice site of ISC of corresponding type. 
It is clear that the all clusters is of the same type, i.e. described by the 
same exponent, which is P/u = 5/48 = 0.10416... as it should be [1] and as 
stressed in [^. A linear fit gives the actual values 0.1056(14), 0.0999(20) and 
0.1046(26), which are not far from the exact value. This is the additional 
argument that we are in force to summarize the probabilities nh^. 

4 Discussion 

We confirm more accurately the simulation result of (g, |^ g for 2D percola- 
tion: the probability of two disjoint Incipient Spanning Clusters has a small 
but finite value. The difference in the values for probabilities is due to the 
fact that in [^, § helical boundary conditions are used. Whereas we used 
free boundary conditions. 

The crossing probability with periodic boundary conditions (PBC) is 



known to be larger and the value 0.63665(8) is computed in ||19[|. We estimate 
in computer simulations the probability of disjoint ISC clusters to be P{k > 
1) « 2.0(4) 10~^ and P{k > 2) ?a 1.4(5) 10"'^ which are even smaller than 
for the case of free boundary conditions (FBC) considered by us. The finite 
size scaling of PBC is more complicated jl^ and fit of the results for only 
four lattice sizes (8,16,32,64) we simulated gives less accurate limiting values 
in comparison with the one for FBC. Surprisingly, the linear fit of logarithm 
of P{k) versus k"^ (shown in Fig. 4 by solid line) gives better accuracy then 
for the case of FBC: apBC ~ 1-915(1) and A ^ 4.26(3). 

It is necessary to note, that we assume the existence of the limiting 
probabilities (2) throughout this work and that the probabilities tt^^ with 
k > 1 reach their maximum in critical region. Indeed, our preliminary 
simulations |2C] shows that the maximum of probabilities tt/u and vr/ij occurs 

^Actually, we got most of the results before getting the last preprint 



at some Pmax, which varies as {pmax —Pc) <x. L^^''^ {v = 4/3 is the correlation 
length exponent) and that the limiting probabilities computed here seems 
to be correct. 

Occasionally, from a total number of 10^*^ samples investigated by us 
(not all the data presented here), we detect the one event of sample with 
just four Incipient Spanning Clusters, which is as probable event as 1/10. 
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Table 1: Probabilities of more than k incipient spanning clusters on critical 
bond square lattices with linear size L and free boundaries. Note, that the 
values of P{k > 1) is multiplyed by a factor 10^ and those of P{k > 2) 
by a factor 10^. For each lattice size L the first row is the probability of 
horizontal crossing and the second one is the probability of vertical crossing. 
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k>0 


fc > 1 10^ 


fc > 2 10^ 


8 


0.50005(5) 


7.657(8) 


3.40(15) 




0.50003(4) 


7.660(8) 


3.98(21) 


12 


0.50002(5) 


7.084(9) 


2.57(14) 




0.49995(5) 


7.070(8) 


2.10(13) 


16 


0.50003(7) 


6.855(9) 


1.97(17) 




0.50002(6) 


6.843(8) 


1.79(19) 


20 


0.49990(6) 


6.742(8) 


1.95(14) 




0.50008(5) 


6.745(9) 


1.72913) 


30 
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1.52(14) 
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1.52(12) 
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0.49999(5) 
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1.73(12) 
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1.56(11) 
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1.51(14) 


oo 


0.50002(2) 


6.58(3) 


1.48(21) 



.. ..Q. .. . 


.. .. Q. .. . 


■ ■ ■ 'O 


... .Q.. .. 


o 


,,,,Q,,,, 


f % 

,,,,Q,,,, 


■ ■ ■ ■© 


... .Q.. .. 


o 


....Q.... 


....Q.... 


■ ■ ■ 'O 


■■■■©■■■■ 


o 


€)■■■■ 


©■ ■ ■ ■ 


■■■■©■■■■ 


..,.Q.... 


....© 


o■■■■ 

(» '= i 


©■ ■ ■ ■ 

► ^ i 


► = i 


..,.0.... 
► < 


t '■ 



Figure 1: Example of simulated lattice with linear size L = 5 (solid lines) 
and its dual (dotted lines). Note, that the number of sites and bonds in 
both directions and for both lattices is just equal to L. 



P{k > 1) • 10^ 
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Figure 2: Probability of more than one Incipient Spanning Cluster multiplied 
by one thousand for 2D bond percolation model as a function of l/L?. The 
linear lattice sizes L are 8,12,16,20,30,32,64. The probability approaches the 
value of 0.00658(53) in the limit of infinite L. Error bars are computed over 
100 bins of 10^ samples each. 
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Figure 3: Probability of more than two Incipient Spanning Cluster multi- 
plied by one million for 2D bond percolation model as a function of 1/L^. 
The linear lattice sizes L are 8,12,16,20,30,32,64. The probability approaches 
the value of 0.00000148(21) in the limit of infinite L. Error bars are com- 
puted over 100 bins of 10^ samples each. 
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Figure 4: Estimated probabilities for more than k Incipient Spanning Clus- 
ters as function of fc^ in logarithmic scale. The dotted line is a fit for lattices 
with free boundaries and the solid line is a fit for the lattices with periodic 
boundaries in vertical direction. 
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Figure 5: Mean size of /c simultaneous Incipient Spanning Clusters, measured 
on the corresponding events. The slope of the curves is close to the exact 
value of fS/u for the all plotted curves. 
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